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ABSTRACT 

We study the effects of internal symmetries on the decay by bubble nucleation 
of a metastable false vacuum. The zero modes about the bounce solution that 
are associated with the breaking of continuous internal symmetries result in an 
enhancement of the tunneling rate into vacua in which some of the symmetries 
of the initial state are spontaneously broken. We develop a general formalism for 
evaluating the effects of these zero modes on the bubble nucleation rate in both 
flat and curved space-times. 



* electronic mail: sasha@langacker.hep.upenn.edu; address after October 1, 1996: Theory Division, 
CERN, CH-1211 Geneva 23, Switzerland 
t electronic mail: klee@phys.columbia.edu 
■'■electronic mail: ejw@phys.columbia.edu 



1 Introduction 



During its early history the universe may have undergone a number of phase transitions. 
One or more of these may have been first-order transitions in which the universe was for a 
time trapped in a metastable "false vacuum" state, with the transition proceeding by the 
nucleation and expansion of bubble of the stable "true vacuum". A crucial quantity for 
the development of such a transition is the bubble nucleation rate per unit volume T. A 
semiclassical procedure, based on a Euclidean "bounce" solution, has been developed for the 
calculation of T at zero temperature H, and extended to the case of bubble nucleation at 
nonzero temperature 0. 

Recently || it was noted that there can be an enhancement of V if a continuous internal 
symmetry of the false vacuum is completely or partially broken by the true vacuum. This 
enhancement can be understood as a consequence of the fact that, instead of a single possible 
final state, there is a continuous family of degenerate true vacuum states into which the decay 
can occur. More formally, the effect arises from the existence of additional zero eigenvalues 
in the spectrum of small fluctuations about the bounce solution. 

The primary focus of Ref. 0] was on the case of a broken U(l) symmetry (see also 0]). 
While a similar enhancement is expected for larger symmetry groups, the treatment of the 
zero modes becomes somewhat more complicated for the case of a non-Abelian symmetry. 
In this note we develop the formalism needed to deal with the case of an arbitrary symmetry. 
We also discuss some further implications of these results, including the extension of these 
results to bubble nucleation in curved space-time. 

The remainder of this paper is organized as follows. In Sec. 2 we develop the general 
formalism and estimate the magnitude of the enhancement that can be achieved. As a 
concrete example, we apply this formalism to the case of SU(2) symmetry in Sec. 3. In 
Sec. 4 we discuss the extension to this work to curved space-time, using the formalism of 
Coleman and De Luccia ||. We show that although the zero mode contribution to the curved 
space-time nucleation rate appears at first sight to be rather different from its flat space- 
time counterpart, it does in fact give the expected result in the the limit where gravitational 
effects are negligible. Section 5 contains some brief concluding remarks. 



2 General Formalism 

We consider a field theory whose fields we assemble into a column vector <fi(x) with purely 
real components. The standard method [|J for the calculation of the quantum mechanical 
(i.e., zero temperature) bubble nucleation rate per unit volume, T, is based on the existence 
of a "bounce" solution <pb{x) of the Euclidean field equations that tends to the false vacuum 
value cf)f at spatial infinity and approaches (although not necessarily reaches) the true vacuum 
value 4> t near a point that may be taken to be the origin. The result may be written in the 
form 

r = -^ (2.i) 

Here If and If, are the contributions to the Euclidean path integral of e~ SE (where Se 
is the Euclidean action) from the homogeneous false vacuum and bounce configurations, 
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respectively, while the division by f2, the volume of the four-dimensional Euclidean space, 
arises in order to obtain a rate per unit volume. The contribution to the path integral from 
a stationary point <j)(x) can be evaluated by expanding the field as 

0(x)=0(x)+]Tc^(x) (2.2) 

3 

where the ipj(x) form a complete set of orthonormal eigenfunctions of the second variation 
operator S'e{4>) = S 2 SE/S(p(x)S(f)(y). For <fi(x) = 0/, this gives, in leading approximation, a 
product of Gaussian integrals over the real variables Cj that results in 

I f = e~ SE ^ {det[S"((j) f )}}- 1/2 = e~ s ^ D f (2.3) 

The calculation of the bounce contribution is complicated by the presence of zero eigen- 
values in the spectrum of S'^(<f)b). Four of these correspond to translation of the bounce 
in the four-dimensional Euclidean space. We will assume that the remainder are all associ- 
ated with internal symmetries of the false vacuum that are not symmetries of the bounce. 
These zero modes are handled by eliminating the corresponding normal mode coefficients q 
in favor of an equal number of collective coordinates Zj. The zero modes about a bounce 
configuration <pb{x, z) can then be written in the form 

Ux,z) = N,,(z)^p^ (2.4) 



where the satisfy the equation 

l{*N)-% = 1^-^ - 2«(M t)u (2.5) 

that follows from the ort honor mality of the ipi. (The factor of 2n is for later convenience.) 
The bounce contribution can then be written as a product of two factors. The first 

j(i) = e -s B (* 6 )I \det'[S"((P b )]\- 1/2 ee e-^W-fc) Db ( 2 .6) 

arises from the integration over the modes with nonzero eigenvalues. The prime indicates 
that the functional determinant is to be taken in the subspace orthogonal to the zero modes, 
while the factor of 1/2 arises from the integration over the single negative eigenvalue mode. 
The second factor, from integrating over the remaining n variables, is 



/f) = (2n)~ n/2 J d n z det 



dcj 
dzj 



(2.7) 



where the factors of 2ir compensate for the absence of n Gaussian integrations. To cal- 
culate the Jacobian determinant, we first equate the change in the field resulting from an 
infinitesimal change in the Zi with that corresponding to a shift of the q, to obtain 



i>i(x,z) dci = d^J' ^ dzj 
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Using the orthonormality of the ipj, we then find that 



(2vr)- n/2 det 

so that 



dcj 
d Zj 



(2tt)"'/ 2 det[M b (z)N^(z)\ = [det M b (z)} 1/2 (2.9) 



if = J d n z[detM b (z)} 1/2 (2.10) 

The fact that the zero modes all arise from symmetries of the theory might lead one to 
expect that the integrand in this equation would be independent of the Zj. Actually, this is 
true only if the measure d n z is invariant under the symmetry transformations. If it is not, 
let fj,(z) be such that fi(z) 1 ^ 2 d n z is an invariant measure and write 

if = J d n zv(z) 1/2 [pi(z)- 1 det M b {z)\ 1 ' 2 (2.11) 

The quantity in brackets is now ^-independent and can be taken outside the integral. 

One can always choose coordinates so that this expression can be written as a product 
of a contribution from the translational zero modes and a contribution from the internal 
symmetry zero modes. For the former, the natural choice of collective coordinates are the 
spatial coordinates z^ of the center of the bounce. The derivative of the field with respect 
to these is, up to a sign, the same as the spatial derivative of the field. Furthermore, with 
these coordinates fi(z) = 1, so the integration over the z^ simply gives a factor of Q. Hence, 
the contribution of these modes to if is f2J^ rans , wheref] 



■h = ( 27r J 



n / d A x(d^f 



1/2 



(2.12) 



The internal symmetry zero modes arise from the action of the gauge group G on the 
bounce solution. Because the bounce tends asymptotically toward the false vacuum </>/, 
normalizable modes are obtained only from the unbroken symmetry group Hf C G of the 
false vacuum. Furthermore, there are no such modes from the subgroup K b C Hf that 
leaves the bounce solution invariant. Hence, the corresponding collective coordinates span 
the coset space Hf/K b . The contribution from these to if is then 

J b Hl/Kb V{H f /K b ) = [ / i(^ )- 1 detM^ / ^(( 7o )] 1 / 2 V{H f /K b ) (2.13) 

where V(Hf/K b ) is the volume of the coset space, M b Hf b is the submatrix corresponding 
to the internal symmetry zero modes, and go is an arbitrary point of Hf/K b . To evaluate 

J b 1 6 it is convenient to take go to correspond to the identity element of Hf. Writing the 
group elements near the identity in the form e iajTj , we may take the collective coordinates 
to be the parameters that multiply the Tj that that span the coset Hf/K b . Evaluated at the 
identity element, the function fi(g) is then equal to unity, while the derivatives with respect 

For the case of a spherically symmetric bounce in a scalar field theory with a standard kinetic energy 
term, J* rans can be expressed in terms of the bounce action, with J* rans = [5b (</>&) — 5e(</)/)] 2 /47t 2 . 
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to the collective coordinates are given simply by the action of the generators on the bounce 
solution. Hence 

1 /2 

J?'' K > = jdet [(27T)" 1 J d 4 x4(x)T]TMx)] } (2-14) 
with the Tj being the generators of Hf/Kb- Gathering our results together, we obtain 

P) Ttrans j^f/^b 

T = e~ B V(H f /K b ) (2.15) 

u s 

where B = S E {4>b) - 5^(0/) . 

It is important to note that Kj, is determined by the symmetry of the bounce, and not by 
that of the true vacuum; it is conceivable (although we believe it unlikely) that the latter is 
invariant under a larger subgroup K t C if/ than the former. Even if K t is identical to K b , it 
is not in general the same as the unbroken symmetry group H t C G of the true vacuum. For 
example, if G is unbroken in the true vacuum and completely broken in the false vacuum, 
Hf, and hence K t , are trivial even though H t = G. In addition, the subgroup K b depends 
not only on the symmetries of the true and false vacua, but also on their relative orientation. 

This last point can be illustrated using a theory with global SU(5) symmetry. Let us 
assume that there is a single scalar field 0, in the adjoint representation, with the potential 
such that the false vacuum has unbroken SU(4) x £7(1) symmetry and the unbroken symmetry 
of the true vacuum is 57/(3) x 57/(2) x U(l). Without loss of generality we may choose the 
false vacuum configuration to be of the form 

(fif = diag (a, a, a, a, — 4a) (2-16) 

The 57/(5) orientation of this field of this configuration influences that of the true vacuum 
bubbles that nucleate within it. Thus, decays to the true vacua with 

0« = diag (MA (2.17) 



and 

# = diag (~b, ~b,b, b, b) (2.18) 

are governed by inequivalent bounce solutions and proceed at different rates ||. If the 
bounces have the maximum possible symmetry, then in the former case K b = ST/" (3) x 
U(l) x U(l) and there are six internal symmetry zero modes, while in the latter K b = 
SU{2) x SU{2) x U(l) x U(l) and there are eight such modes. Of course, bubble nucleation 
could with equal probability lead to any configuration obtained by applying an SU (4) x U(l) 
transformation to (ffp or ; this is taken into account by the integration over the coset 
spaces (£17(4) x U(1))/(SU(3) x 17(1) x 17(1)) and (517(4) x U{l))/{SU{2) x 517(2) x *7(1) x 
U(l)). However, there are true vacuum configurations that cannot be obtained by such 
transformations. In general, there are no bounce solutions corresponding to these, reflecting 
the fact that if a bubble of such a vacuum were to form, the external false vacuum would 
exert forces that would realign the field in the bubble interior. 

Finally, let us estimate the magnitude of the zero mode corrections that we have found. 
For defmiteness, we will consider the case of a scalar field theory whose potential can be 
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written as V{4>) = AF(0) with F((f>) containing no small dimensionless parameters. Standard 
scaling arguments using the fact that the bounce is a stationary point of the action show 
that the bounce has a radius ~ 1/m (where m is a characteristic scalar mass) and an action 
(relative to that of the false vacuum) of order 1/A. The typical magnitude of the bounce field 

is (fib(x) ~ m/y/X, while the Tj are all of order unity, so J b Hf ^ Kb ~ (Am 2 ) - ^ 2 , where N is 
the number of internal symmetry zero modes. The ratio D b /Df is of order unity in A, but is 
proportional to a dimensionful parameter ~ m N+A arising from the fact that the contribution 
of the zero eigenvalue modes has been deleted from Df. Finally, the coset volume is of order 
unity. Overall, then, we have 

r = Cl A-( w '/ 2 mV C2 / A (2.19) 

where q and c 2 are of order unity; the effect of the internal symmetry zero modes has been 
to enhance the nucleation rate by a factor of order \~ N I 2 . Phrased somewhat differently, the 
enhancement is roughly by a factor of B N I 2 . 



3 SU{2) Symmetry 

As a concrete example, let us consider the case where the symmetry group of the false 
vacuum is Hf = SU{2) but the bounce solutions break this symmetry. A natural set of 
collective coordinates is given by the Euler angles. Thus, given one bounce solution 4>l(x), 
we can define a three-parameter family of solutions by 

<f> b (x; V, 0, = <fi(x) = U(<p, 9, ^)<f>° b (x) (3.1) 

where the Tj are the appropriate (possibly reducible) representation of the generators of 
577(2). Differentiation of this expression gives 

d v <f> b (x;<p,9,i>) = iU(<p,9,i>)f 3 <f> b (x) 

fy&fo^M) = iU(ip,6^)T3<f>° b (x) (3.2) 

where 

= cos ip sin 9 Ti + sin ip sin 9 T 2 + cos 9 T 3 
T 2 = e _ *^ T3 T 2 

= — sin ijj Ti + cos ip T 2 (3.3) 

Thus, if Zj = ((p,9,tp), 

d j( f> b (x; tp, 9, i>) = iK jk U(<p, 9, ^)T k 4>l{x) (3.4) 

where 



cos ij) sin 9 sin tp sin 9 cos ip cos 9^ 
K = | — sin^ cos^ 

I J 



(3.5) 
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Substitution of this into Eq. (|2.5|) yields 

{M sum/K h)u = {27t) -i k . jKh J <? x $( x )T}T k <t>l{x) (3.6) 

Now recall that an invariant measure on SU(2) is given by sin 8d(f d9 dip, so we may take 
/i(y?, 8, ip) = sin 2 8 = (det Kf . Hence, 



fi- 1 det M^ u{2)/Kb = det (27T)- 1 J d A xcjP b \x)T}T k (jP b {x) 



(3.7) 



This is independent of the collective coordinates, as promised, and is in agreement with 
Eq. (gig) . 

Three specific cases may serve to illustrate some of the possible behaviors: 

a) One SU(2) doublet: If the bounce involves an 577(2) doublet, then the bounce com- 
pletely breaks the SU (2) symmetry. The coset volume factor is 

V(H f /K b ) = V(SU{2)) = 16tt 2 (3.8) 
while Mfr U b = M b U is proportional to the unit matrix, with 

= < J ii (27r)- 1 / d 4 x4(x)<j) b (x) (3.9) 



A*" 1 



M SU(2) 



(Because our formulas have been derived using real fields, one must use a four-dimensional 
real representation, rather than a two-dimensional complex representation, for the Tj when 
obtaining this result from Eq. ( |3.7| ).) 

b) One SU(2) triplet: If the bounce is constructed from a single real triplet whose 
direction is independent of x [i.e., such that (j) b {x) can be written as (0, 0, f(x))}, then K b = 
£7(1). There are only two zero modes and Hf/K^ is the two-sphere spanned by 8 and tp with 

V(H f /K b ) = V(SU(2)/U(1)) = 4vr (3.10) 

j\^su{2)/K b . g a 2 x 2 matrix, with 



M SU(2)/U(l) 



>J 



S^n)- 1 J d A xct> 2 b {x) (3.11) 



c) Two non-parallel SU{2) triplets: If the bounce solution contains two triplet fields that 
are not parallel, then the bounce has no continuous symmetry. Because only integer spin 
fields are involved, 

V{H f /K b ) = V(SO(3)) = 8tt 2 (3.12) 
The matrix n^M^ ^ has three unequal eigenvalues. 

4 Bubble Nucleation in Curved Space-Time 

Coleman and De Luccia 0] showed that the bounce formalism could be extended to include 
the effects of gravity by requiring that both the bounce and the homogeneous false vacuum 
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configurations be solutions of the coupled Euclidean matter and Einstein equations. For 
a scalar field theory with V(4>) > 0, as we henceforth assume, the false vacuum solution 
consists of a uniform scalar field ^ on a four-sphere of radius 



3M « (4.1) 



with total Euclidean action (including gravitational contributions) 

s ^ = -wm < 4 ' 2 > 

The bounce solution has the same topology, with regions of approximate true vacuum and 
false vacuum separated by a wall region. If the matter mass scale M. is much less than 
the Planck mass Mpi, then both the radius Rf, of the true vacuum region and the difference 
between the bounce action and the false vacuum action differ from the corresponding flat 
space quantities by terms of order (.M/Mpi) 2 . 

The spectra of the small fluctuations about these solutions again contain one zero mode 
for each symmetry of the Lagrangian that is broken by the solution. However, because the 
Euclidean solutions are on closed manifolds with finite volumes, the modes due to symmetries 
broken by the false vacuum are normalizable, in contrast with the flat space case. Hence, 
we would expect the flat space factors given in Eq. ( 2.13Q to be replaced by 



J? ,Kb V{G/K b ) _J^ 

Although the volume factors give the same result as in flat space, the Jacobean factors 
appear quite different. Yet, for M. <C M P1 , where gravitational corrections should be small, 
this should approach the flat space result. 

To see how this comes about, let us denote by tj the generators of Hf/Kj, and by Sj those 
of G/Hf. The Jacobean determinant in the numerator of Eq ( |4.3| ) has contributions from 
matrix elements containing both types of generators, whereas the determinant in the denom- 
inator only involves SiSj matrix elements. Because the tj annihilate the false vacuum, the 
matrix elements involving these have nonzero contributions only from the region, of volume 
~ Rf, where the bounce solution differs from the false vacuum and hence are suppressed by a 
factor of order (HfRb) 4 ~ (A^/Mpi) 4 relative to the s^Sj matrix elements. (We are assuming 
that Rf, ~ Ai^ 1 ; this will be the case for generic values of the parameters.) This implies 
that, up to corrections of order (J\4/Mp{) 8 , the determinant can be written as a product of 
a determinant involving only the tj and one involving only the Sf, i.e., 

jG/K b jG/Hf 

im = -7m J " f/Kb t 1 + °((^/ m pi) 8 )] ( 4 - 4 ) 

The first factor on the right hand side differs from unity only by an amount proportional 
to the fraction ~ (.M/Mpi) 4 of the Euclidean space where the bounce differs from the 
false vacuum. The second factor differs from the corresponding flat-space term only by the 
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replacement of the matter fields of the flat space bounce by those of the curved space bounce, 
and so clearly reduces to the flat space result as Ai/M P \ — > 0. 

The fact that the bounce solution is a closed manifold, with the true and false vacuum 
regions both finite, suggests that it can contribute not only to the nucleation of a true vacuum 
bubble within a false vacuum background, but also to the nucleation of a false vacuum bubble 
within a true vacuum background, with the rate for the latter process obtained from that of 
the former by making the substitution 0/ — > <p t [/J. To leading order, the ratio of these two 
rates is 



iS B ((l>t)-S E (<t>f) 



exp 



1 



(4.5) 



The continued nucleation and expansion of bubbles of one vacuum within the other will 
result in a spacetime that is a rather inhomogeneous mixture of the two vacua. There is an 
intriguing thermal interpretation of this mixture if V(<pf) — V(<p t ) <§C (V (<p /) +V (<fit)) / '2 = V, 
so that the geometry of space is approximately the same in the regions of either vacua, with 
a Hubble parameter 



H 



8ttV 



\ 3M] 



(4.6) 



It seems plausible that the fraction of space contained in each of the vacua might tend to 
a constant, with the nucleation of true vacuum bubbles in false vacuum regions being just 
balanced by the nucleation of false vacuum bubbles in true vacuum regions. For such an 
equilibrium to hold, the volumes Qf and Q t of false and true vacuum must satisfy 



5/ 
fit 



r 



-n boI [v(<t> f )-v(<t>t)]/T H 



(4.7) 



where the horizon volume i\or = (47r/3)i/~ 3 and the Hawking temperature T# = H/2n. If 
we view the de Sitter space as being somewhat analogous to an ensemble of quasi-independent 
horizon volumes in a thermal bath, then this leading contribution to the volume ratio is 
essentially a Boltzmann factor. 

The zero mode corrections to the nucleation rate are consistent with this thermodynamic 
picture. Their effect is to multiply the ratio in Eq. ( |4.7|) by 



V(G/H f ) J[ 



G/Hf 



H 



N t -N f 



V(G/H t ) jO/B t 



7T 



V{G/H f ) 
V(G/H t ) 



det 



det 



(JWZh/27t) (4>)TiT r 



(n hoi T H /27r) UlTiTj 



1/2 



(4- 



where Nf and N t are the number of internal symmetry zero modes in the false and true 
vacua, respectively. We recognize the dimensionless ratio on the right hand side as the ratio 
of two classical partition functions of the form 



d N zd N p _ 



H Z /T H 



d N z 



{2ti) n 

that follow from the effective Lagrangian 



(n hor T H /2n) (pTjTj 



1/2 



J'T I Z i Z j 



(4.9) 



(4.10) 
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that describes the collective coordinates dynamics for a horizon volume with spatially uni- 
form scalar field <j). 

The presence of a dimensionful prefactor in Eq. (|4.8|) is required by the differing numbers 
of zero modes about the true and false vacua, which implies a dimensional mismatch between 
the functional determinants over the nonzero eigenvalue modes. This suggests that the factor 
of (H / \/3 %) Nt ~ N f should, like the functional determinants themselves, be understood as 
related to the first quantum corrections to the vacuum energies. 

5 Concluding Remarks 

We have seen that when a metastable false vacuum decays to a true vacuum that breaks 
some of the internal symmetries of the false vacuum, the presence of iV > zero modes 
about the bounce solution can lead to an enhancement of the bubble nucleation rate. In 
a theory characterized by an overall scalar coupling A, the zero-temperature, quantum me- 
chanical, tunneling rate is increased by a factor of order \~ N ^ 2 ~ B N I 2 . A straightforward 
extension of our methods to finite temperature thermal tunneling shows that, although the 
A-dependence is changed, the enhancement is still of order B N I 2 . Since the nucleation rate 
falls exponentially with B, we therefore have the curious situation that the enhancement is 
greatest when the overall rate is smallest. 

These results may be of particular interest for the symmetry-breaking phase transitions 
that arise in the context of a grand unified theory. For any given nucleation rate, the numer- 
ical effect of the zero mode corrections will, of course, almost always be negligible compared 
to the uncertainties due to the undetermined couplings in the scalar field potential. The zero 
mode effects could, however, be significant when there are competing decays to vacua with 
different degrees of symmetry breaking, such as are encountered in many supersymmetric 
models. 
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